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Characterizing the satellites of massive galaxies up to z ~ 2: young 
populations to build the outskirts of nearby massive galaxies 
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ABSTRACT 

The accretion of minor satellites is currently proposed as the most likely mechanism to 
explain the significant size evolution of the massive galaxies during the last ~ 10 Gyr. In 
this paper we investigate the rest-frame colors and the average stellar ages of satellites found 
around massive galaxies (M star ~ lO n M ) since z ~ 2. We find that the satellites have bluer 
colors than their central galaxies. When exploring the stellar ages of the galaxies, we find that 
the satellites have similar ages to the massive galaxies that host them at high redshifts, while 
at lower redshifts they are, on average, > 1 .5 Gyr younger. If our satellite galaxies create the 
envelope of nearby massive galaxies, our results would be compatible with the idea that the 
outskirts of those galaxies are slightly younger, metal-poorer and with lower [a/Fe] abundance 
ratios than their inner regions. 
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1 INTRODUCTION 

There is a growing consensus that the str ong size evolution 
observed among the mas s ive galaxies (e.g.. iDaddi et al.l 120051: 



Truiill o et alj l2006t 120071 : iLonghetti et al.l 120071 : iBuitrago et al.l 



2008) is mainly dominated by the continuous accretion of satel 



lites, which would increase the size (by a factor of ~ 4) and the 
mass (by a factor of ~ 2) of the galaxies since z ~ 2. The theoret- 
ical expectations associated to this mechanism sugg ests a growth 



in size compatible with the observations (e.g. iNipoti et al 



Khoc hfar & Burkertll20od : iHopkins et alll2009l : iNaab et al 



2003 



2009 



Oser et al.ll2012l) . In this scheme, the infalling satellites would rep- 
resent the building blocks of the external parts of the most massive 
galaxies. 

A direct way of probing the above minor merging scenario 
is to explore the evolution o f the satellites around massive galax- 
ies with cosmic time (e.g.. IWilliams et al.l 1201 ll : iNewman et al.l 
l2012l : IOuilis & Traiilldl2012h . Several papers in the last few years 
have calculated the frequency of satellites around massive galaxies 
and have quantified how this fraction changes with cosmic time. 
This work has been d one in the nearby Universe (see e.g., Ichenl 
2008: lLiuet"aT]|201lh and up to z ~ 1 usi ng mostly samples of 
central galaxies less m assive than 10"M G (de Ravel et al] |201ll : 
iNierenberg et al.l 1201 ll) . Recently, iMarmol-Oueralto et al.N2012h 
(hereafter Paper I) expanded the previous analysis exploring the 
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fraction of massive galaxies (M stm ~ 10" M Q ) with satellites up 
to z ~ 2. They probed two different redshift ranges: up to z = 1, 
to explore the fraction of massive galaxies with satellites having 
0.01 < M sat /M ccntral < 1 (1:100), and up to z = 2, to study the 
fraction of massive galaxies with satellites within the mass range 
0.1 < M sat /M ccntra i < 1 (1:10). The results of that paper support 
a picture where the fraction of massive galaxies with satellites, 
within a projected radius of 100 kpc, has not changed with time 
since z ~ 2. This fraction remains around ~ 15 % for galaxies 
with satellites with mass M star > 10 10 M G and around ~ 30% if we 
explore satellites with masses M star > 10 9 M o up to z = 1. 

In this paper we investigate the stellar population properties 
of the satellites around massive galaxies found in Paper I in an ef- 
fort to characterize the satellites that eventually could be accreted 
to form the outskirts of the massive galaxi es that we find in the 
local univers e . Using the Rainbo w database JPerez-Gonzalez et al.l 
120051 12008 j: iBarro et al 1 l201 lah we analyze the evolution of the 
rest-frame colors, as well as the average stellar ages of the satel- 
lites in comparison with their central galaxies. 

This paper is structured as follows. In Section [2] we describe 
our sample of massive galaxies and the photometric catalog we 
have used to identify their satellites and their properties. A sum- 
mary of our criteria to identify the satellites is presented in Sec- 
tion l2.ll We present in Section[3]the results of this work related to 
the rest-frame colors l[3}, color-magnitude diagrams 13.2b and av- 
erage stellar ages d3 - 3t . Finally, a summary and discussion of our 
findings is provided in Section [4] In this paper we adopt a stan- 
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dard ACDM cosmology, with _2 m = 0.3, fi A = 0.7 and H = 70 
km s _I Mpc~'. All magnitudes are in the AB system (Oke 1974). 



2 THE DATA 

We have used as the reference catalog for the central galaxies the 
compilation of 629 massive (M stal . ~ 10" M Q ) galaxies presented 
in Paper I. These galaxi es were selected fro m the sample of mas- 
sive objects published in lTruiillo et alj J2007l) (h ereafter T07) based 
on the near-infrared Palomar/DEEP-2 survey teundv et al1l2006l : 
IConselice et al.l2007h over 710 arcmin 2 in the Extended Groth Strip 
(EGS). T0 7 estimated ( circularized) half-light radius (r e ) and Sersic 
indices n ( Sersic] I"l968l) using the ACS /-band filter for all the 
galaxies in our sample. These authors tested the accuracy of the 
parameter determination running a set of simulations of artificial 
galaxies uniformly generated at random, matching the observed 
structural parameter distribution of their galaxies. Moreover, they 
simulated the real conditions of their observations (the background 
sky and the observed PSF). They obtained that the uncertainties of 
the structural parameters were Sr c /r c < 30 per cent and Snjn < 38 
per cent. 

The satellites around our massive objects were extracted from 
the EGS IRAC-selected galaxy sa mple from the Rain bow Cos- 
mological Databasjj] published b y iBarro et al] fcOllal) (see also 
IPerez-Gonzalez et al.l2 005. 2008a). This database covers an area of 
1728 arcmin 2 centered on the Extended Groth Strip (EGS) and pro- 
vides spectral energy distributions (SEDs) ranging from the UV to 
the MIR regime plus we ll-calibrated and re liable photometric red- 
shifts and stellar masses tearroetalj2011bl) . For the ~ 10% of the 
galaxies in this catalog, spectroscopic redshifts are also available. 
From the Rainbow database we have selected all the galaxies with 
z < 2.2 and an estimated stellar mass 10 8 M Q < M 8tar < 10 I2 M o . 
A total of ~55000 objects were selected in the whole EGS area 
following these criteria. We refer to this resulting sample as the 
Rainbow catalog. 

There are 629 galaxies up to z = 2 from T07 for which the 
Rainbow catalogue allows the study of satellites down to a 1:10 
mass ratio (0.10 < M sat /M centra i < 1). Down to 1:100 mass ra- 
tio (0.01 < M sat /M ccntra i < 1) and z < 1, the number of galax- 
ies that can be explored is 194. These final samples are given by 
the stellar mass limit (75% complete) at each redshift of the Rain- 
bow database (see Fig. 1 of Paper I), so that we only consider cen- 
tral galaxies that could have satellites down to 1:10 (and 1:100 for 
z < 1) mass ratio in the Rainbow catalogue. More information of 
this selection is detailed in Paper I. Briefly, we compared the data 
for the massive galaxies (redshifts and stellar masses) from T07 
and the Rainbow catalog to build a sample of 694 massive galax- 
ies, of which 317 have photometric redshifts from Rainbow and 
377 have spectroscopic redshifts (340 from the Rainbow catalog 
and 37 from T07). After this selection, we applied a final cut in this 
sample of galaxies to assure that the fraction of galaxies with satel- 
lites along our explored redshift range is not biased by the stellar 
mass completeness limit of the Rainbow database (see their Fig. 4 
Perez-Gonzale z et al. 2008a) in the redshift range < z < 2. The 
mean stellar mass of our sample is M 4 = 1.3 x 10" M Q (Kroupa 
IMF). 

Finally, our sample of massive galaxies is large enough that 
we can explore whether the properties of the satellites around the 
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massive galaxies depends on the morphological type of our ob- 
jects. We have used the Sersic index as a proxy to the galaxy 
morphology, since in the nearby universe, galaxies with n < 2.5 
are mostly disc-like o bjects, whereas galaxies with n > 2.5 are 
mainly spheroids (e.g., Andredakis et al.ll995l : lBlanton et alj2003l : 



iRavindranath et alj|2004h . We have used the published Sersic in- 
dexes provided by T07 separate our galaxies in two sub-samples: 
spheroid- and disk-like galaxies, in our whole redshift range. 

2.1 Search for satellites around the massive galaxies 

In this work we will use the results already presented in the Table 1 
of Paper I related to the fraction of massive galaxies with satellites. 
Summarizing, we counted as satellites those galaxies found in the 
Rainbow catalog that: (1) are within a projected radial distance to 
our central galaxies of i? seilr ch = 100 kpc (corresponding to 0.3 and 
0.2 arcmin for z = 0.5 and 2.0, respectively). The search is also 
restricted to distances larger than 1 arcsec (~ 8 kpc), which is the 
deblending limit of sources in the Rainbow database; (2) the dif- 
ference between their photometric redshifts and the redshift of the 
central galaxies is lower than the lcr uncertainty in the estimate of 
the photometric redshifts of the Rainbow database (using the same 
uncertainties than in the selection of the sample, i.e., Az p hot = 0.070 
for 0.0 < z < 0.5, Az phol = 0.061 for 0.5 < z < 1.0, and 
Azphot = 0.083 for 1.0 < z < 2.5); and (3) the stellar mass of these 
objects should be within 0.1 < AZ sat /M centla i < 1.0 for the galaxies 
in the range < z < 2, and within 0.01 < M sat /M cclltral < 1.0 for 
the galaxies in the range < z < 1. We considered different red- 
shift bins to explore the evolution of the fraction, F sat , of massive 
galaxies with satellites. The width of these bins (see Table[T} were 
chosen to include a similar number of massive galaxies in each bin 
and have a similar statistics among them. 

Although redshifts, either from photometric or spectroscopic 
measurements, are available for all the galaxies in this study, there 
is a number of objects identified as satellites that are actually con- 
taminants that satisfy the above criteria but are not gravitationally 
bound to our massive galaxies. It is important to note that in ulti- 
mate case it is not the accuracy of the redshift estimate which lim- 
its the determination of pure gravitationally linked system s (e.g., 
IPatton et alj200d ; |Patton & Atfieldll2008l ; lLotz et all201ll) . but the 
deviation from the Hubble flow associated to the peculiar velocities 
of the galaxies. This effect is particularly relevant in dense enviro- 
ments as clust ers of galaxies , and its constribution is difficult to 
estimate (e.g., iLiu et al]|201ll . Paper I). 

To quantify the contamination, we follow the method pre- 
sented in Paper I. Briefly, the method consists on placing a number 
of mock massive galaxies (equal to our central galaxies) randomly 
through the volume of the catalogue. In our simulations, the num- 
ber of mock galaxies that are within each redshift bin is the same 
than in our sample and we keep fixed the parameters of the massive 
galaxies. Once we have placed our mock galaxies through the cata- 
log, we count which fraction of these mock galaxies have satellites 
around them taking into account the searching criteria explained 
above. This procedure is repeated one million times to have a ro- 
bust estimate of the fraction of mock galaxies with satellites. We 
consider this value, 5 s i mu i, to be representative of the background 
affecting our real central sample. In addition, these simulations al- 
lowed us to estimate the scatter in the fraction of galaxies that have 
contaminants. We use this scatter as an estimate of the error of our 
real measurements. 

Finally, the fraction of galaxies with real satellites were com- 
puted following the expression (see Paper I): 
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(1) 



where F sat is the final (corrected) fraction of galaxies with satellites 
obtained from the observed values F obs and the fraction of galaxies 
with satellites obtained in the simulations 5 s i mu i. 



2.2 Rest-frame colors 

For all the galaxies in the Rainbow catalogue, including the massive 
galaxies, there are available measurements of their their rest-frame 
colors. These quantities were obtained from the fitting to the ob- 
served UV-to-MIR photometry (from the ultraviolet wavelengths 
probed by GALEX to the far-infrared observed by Spitzer) com- 
piled for each galaxy in the Rainbow catalogue. The fit is based on 
a set of empirical templates comp uted from PEGASE 2.0 m odels 
dFioc & Rocca-Volmerangell 19971) assuming a lSalpeteil d 19551) ini- 



tial m ass function (IMF; M stiu = 0.1 - 100M m ). a ICalzetti et ail 
(2000) extinction law, and the FIR photometry (M I PS 24 and 
70 iim) of dust emi ssion models of (chary & Elbazl {200 lh and 
bale & Helov] d2002h (we refer the reader to lBarro et aU201 IB . for 

specific details of the se fittings). 

With these data, iBarro et ail J201 lbl) provide with rest-frame 
colors in different filters, not corrected by internal extinction and 
measured over the best-fitted templates of the SED for each galaxy. 
Here we will use the colors u',g', and r 1 of the SDSS system, and 
K s of 2MASS, given by the Rainbow database for each galaxy. 



3 RESULTS 

To estimate properly the average properties of the satellite galaxies, 
we need to make statistical corrections taking into account that a 
substantial part of our satellites are background contaminants. As 
explained in Paper I, this can be done by using the expression: 

^ ^obs ^simul s-. 
< 2sat >= < Gobs > ~ < Gstaul > (2) 

r sat **sat 
where < gobs > is the observed mean value of the property Q (in 
this work, rest-frame colors and stellar ages), < gsimui > is the mean 
value obtained from the simulations (i.e. the values that are found 
for the contaminants) and < Q sat > is the value after the correction. 

For the central galaxies with satellites, it is also necessary to 
take into account that it is a priory unknown which central galaxies 
have real satellites and which ones have only fake satellites. Since 
this is the same problem that we face for the properties of the satel- 
lites, again we will use the Eq. [2] to estimate the mean properties 
of the central objects. In this case, to model the mean properties of 
the galaxies with fake satellites we use the mean properties of the 
central galaxies without satellites. For both the satellites and their 
central galaxies, two sources of error are considered: the error in 
their mean values, and the error in the estimate of the fraction of 
galaxies with satellites. 

Finally, for the massive galaxies without detected satellites, 
we simply compute the mean value of their properties for all the 
galaxies in each redshift bin. The associated errors in this case are 
the errors in their mean values. 



3.1 Evolution of the rest-frame colors with redshift 

We study in this section the evolution of the average colors of the 
central galaxies and their satellites. We focus on u' - r', which 
mainly traces young stellar populations, and g' - K, which takes 



into account the evolved stellar populations. The results are shown 
in Fig.Q] 

We find that massive galaxies with and without satellites have 
similar u' - r' and g' — K colors, and they are redder than their 
satellites. This is in agreement with previous works in the local 
universe (e.g., IWeinmannetal]|2009l using data from the Sloan 
Digital Sky Survey, also pointed to bluer colors for satellites than 
central galaxies in our stellar mass ranges) and up to z = 1 (e.g., 
iGrutzbauch et al.l201 lL using DEEP2 data). We see that at low red- 
shift, our satellites have redder u' - r' (and g' - K) colors than the 
field g alaxies, in agreement with recent works in the loc al universe 
(e.g.. Ivan Den Bosch et ai]|2008l ; IWang & White! |2012|) . while at 
higher redshifts our satellites and the field galaxies present similar 
colors. 

When we segregate the central galaxies depending on their 
morphological type, we see that the massive spheroid-like galax- 
ies have redder u' - r' colors than their satellites, while they present 
similar g' - K. Moreover, disk-like massive galaxies and their satel- 
lites, in general, are bluer that the massive spheroid-like galaxies 
and their satellites. 

In these figures we have also included the expected col- 
ors for the passive evolution of the galaxies. For doing that, we 
measure the color s over the single stellar population models of 
ICharlot & Bruzuall d2009) at different ages, assuming solar metal- 
licity and a Kroupa IMF, and we chose that model with closer color 
to our galaxies at higher redshift with satellites. Then, we plot the 
color evolution from that initial value. In addition, the errors in the 
colors were also considered, giving the upper and lower limits rep- 
resented in the figure. If we compare our data with these colors, we 
find that the u' - r' of the central galaxies at lower redshifts are not 
compatible with a passive evolution. This is not surprisingly, since 
the u' - r' color is mostly tracing the young stellar population that 
could be formed along their life time. However, the passive evolu- 
tion of the g' — K values is more compatible with the colors of the 
central galaxies since this color is showing the evolved stellar pop- 
ulation that was born at similar ages. The evolution of - r 1 and 
g' — K colors point to that new star formation is been produced at all 
redshifts among the different galaxy types, includ ed the spheroid- 
like objects (see also lPerez-Gonzalez et al . 2008^)- 



3.2 Color-magnitude diagrams 

In the local Universe, it is well known the bimodality of galax- 
ies found in the color-magnitude diagram, where early-types 
(E/S0) populate a narrow red sequence th at is separated from 

11 20011; 



bluer, star-forming spirals by a green v a lley {Strateva et al.l 200 ll 
iHogg et alj l2003t iBalogh et alj |2004 iBaldrv et al.l 12004 and 
references t herein). A sim i lar division extends at least up to 
z ~ 1 (e.g.lLinetalJ|l999l;lBelletalJl2004l ; IWeineretalJl2005l; 
Will mer et alj2006l:IWang et al.l 2010) and beyond dGiallongo et ar 
20051 : iPannella et al.l 120091 ; iBrammer et alj 120091 IWhitaker et al 
l201lh . Here we present the color-magnitude diagram for the mas- 
sive spheroid- and disk-like galaxies in our sample in comparison 
with their satellites to explore whether our objects are in the red 
sequence or in the blue cloud and to probe whether there is hint of 
any evolution. 

In Fig[2]we present the u' - r' rest-frame color vs. M r dia- 
gram in the three redshift ranges of study. Since the separation be- 
twee n the red seque nce and the blue cloud evolves with redshift 
(e.g. lBelletalJl2004h . we also present the histograms of the galax- 
ies in the Rainbow database in each redshift range to identity the red 
sequence and the blue cloud in our data. We fit a double Gaussian to 
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Figure 1. Evolution of the rest-frame u' - r' and g' — K colors with redshift. The mean colors of central galaxies with (magenta squares) and without (cyan 
triangles) satellites are plotted together with the mean colors of the detected satellites (blue circles) and the field galaxies (gray stars). Gray shadows represent 
the passive evolution based on the colors of the average population of the massive galaxies with satellites at the higher redshift. 



these histograms as a first order approximation to obtain these two 
populations in our sample. This separation is showed as a dashed 
line in Fig(2] We find that for galaxies with z < 1.10, the mas- 
sive central spheroids and their satellites are in the red sequence, 
although there is a trend of satellites getting bluer as redshift in- 
creases. Also the massive disk-like galaxies at z < 1.10 are mainly 
in the red sequence, while their satellites are in the blue cloud. Fi- 
nally, the satellites of both the spheroid-like galaxies and the disks 
at 1.10 < z < 2.00 are in the blue cloud, while the massive galax- 
ies, independently of their type, are in the red sequence. For disk 
galaxies the results are very uncertain due to the low statistics. 

These color-magnitude diagrams suggest that mergers involv- 
ing massive objects at higher redshifts (1.10 < z < 2.00), would 
be mixed mergers (i.e., a massive galaxy in the red sequence and a 
satellite in the blue cloud), for both spheroid- and disk-like central 
galaxies. Also we could expect mixed mergers for disk-like galax- 
ies at lower redshifts. Finally, the satellites of massive spheroids at 
0.2 < z < 0.75 are more located towards the red sequence, and 
therefore, mostly dry mergers would be expected at those redshifts. 

3.3 Stellar ages of the satellite galaxies 

We investigate in this section the average stellar ages of both mas- 
sive and satellite galaxies. We take the average stellar ages from the 
best fitting templates to the SEDs of each galaxy computed as part 
of the Rainbow database. For this experiment, the templates were 
computed from a pure exponential-declining star formation history 
dSandadfl986l) : 

SFR (t) = exp(-(r - 0/t) (3) 

where the steepness of the decay are regulated by a single param- 
eter t. Here we considered a simple (not delayed) exponential, be- 
ing T = 13.7 Gyr the age of the stellar population of the galaxy 
today. We considered s ingle stellar popula t ions w ith ages ranging 
10 6 - 2 x 10 10 yr from lCharlot & Bruzuall j2009h and t = 8 - 12, 



Table 1. Stellar age differences between central and satellite galaxies. 



Redshift range 


Agecentral " A g e satcllite 




(Gyr) 


All galaxies 



0.10 < M sat /M centra , < 1.00 
0.20 < z < 0.75 0.92 ± 0.06 

0.75 <z< 0.90 1.04 ±0.13 

0.90 <z< 1.10 0.37 ±0.07 

1.10 <z< 2.00 -0.20 ±0.15 



0.01 < M sat /M centra i < 1.00 
0.20<z< 0.55 1.76 ±0.33 

0.55 <z< 0.73 1.11 ±0.07 

0.73 <z< 1.10 0.83 ±0.06 



Spheroid - like (n > 2.5) galaxies 

0.10<M sat /M centra i < 1.00 
0.20 < z < 0.75 0.99 ± 0.06 

0.75 <z< 1.10 1.00 ±0.04 

1.10 <z< 2.00 0.01 ±0.15 



Disk - like (n < 2.5) galaxies 



0.10 < M sat /M ccntra i < 1.00 
0.20 < z < 0.75 0.41 ± 0.43 

0.75 <z< 1.10 -0.73 ±0.30 

1.10 <z< 2.00 -0.19 ±0.42 



and we assumed a iKroupal feOOlh IMF and solar metallicity. A 
dust attenuation ranging A v = 0.0 - 2.0 mag is allowed to fit 
the SED of each galaxy. Here, we analyze the average stellar ages 
of the templates that better fits the SEDs. It is worth noting that 
photometrically-derived ages are not very robust due to degenera- 
cies with tau and metallicity and for that reason we do not intend 
here to provide absolute determination for the ages of the galaxies 
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Figure 2. Color-magnitude diagram for the galaxies in this study, showing 
u' - r' rest-frame colors vs. r-band absolute magnitudes. Red circles rep- 
resents massive spheroid-like galaxies, while blue stars represents massive 
disks. Their satellites are represented with squares of the same colors. The 
dashed magenta line corresponds to the separation for the red sequence and 
the blue cloud in our sample, obtained from a double Gaussian fit to the 
histograms of the Rainbow data in each redshift range, presented on the 
right. 



in this study, just to have a hint of their relative values between the 
central galaxies and their satellites. 

We summarize in Table [T] and Fig. [3] the difference in age of 
the central galaxies and their satellites for the redshift bins consid- 
ered in this study. In all the cases, we find that galaxies at higher 
redshifts present similar average ages, while satellites are younger 
than their central galaxies when the redshift decreases. This ten- 
dency is clear for galaxies with satellites in both 1:10 and 1:100 



mass ratios (left panel in Fig.O. A linear fit indicates that satellites 
at z = would be > 1.5 Gyr younger than their central galaxies. 

When we separate the sample depending on their morpholog- 
ical type (right panel of Fig. [3}, we find the same result for the 
satellites around spheroid-like galaxies, being > 1.5 Gyr younger 
at z = 0. This is not the case for disk-like galaxies and their satel- 
lites which present similar ages within the errors at all the redshifts. 



4 SUMMARY AND DISCUSSION 

Simulations of minor merging suggest that new accreted stars from 
satellites are mainly added to the outer part s of the central galax - 
ies (e.g jNipoti et al.ll2003l : lNaab et all2009l ; lHopkins et alj201Ch . 
This pr oduces an efficient grow in size compatible with the obser - 
vations faopkins et al.l2009l ; lNaab et alj2009l ; lTruiillo et alj201ll) . 
If the satellites are significantly different in stellar content than their 
central galaxies, the imprints of these satellites might be found in 
the outer parts of massive galaxies in the local universe. In this 
paper we explore several properties of the satellite galaxies in com- 
parison with their central galaxies to characterize the galaxies that 
eventually could be accreted by the massive galaxies at different 
redshifts. 

Our sample of central galaxies is formed by 629 massive 
(M sar ~ 10"Mq) ga l axies at 0.2 < z < 2.0 presented in 
Marmol-Oueralt o et al.l d2012h . We searched for satellites within 
a projected radial distance of 100 kpc around the massive galax- 
ies that have mass ratios of galaxies in two regimes: 0.01 < 
M sat /Af central < 1 (1:100) up to z = 1 and 0.1 < M sat /M centra , < 1 
(1:10) up to z= 2. 

When exploring rest-frame colors, we find that satellite galax- 
ies are bluer than their central galaxies, with and witho ut satellites. 
This r esult is in agreement with the recen t works by Wang et al. 
j2010h (data from the G OODS survey iGiavalisco et alj l2004h 
and Newman et al. | j2012h (results from the CANDELS survey; 
Groni n et alj201 1 ), who also explored the colors of satellite galax- 
ies around galaxies. Although with different samples and analy- 
sis, their studies agree in the bluer colors for the satellites when 
comparing with their central galaxies. Furthermore, the satellites 
at low redshift have redder u' - r' colors than galaxies in the field 
(Fig. IB • Several works in the local universe also point out to redder 
colours for the satellites than for galaxies in the field. For example, 
using data from th e Galaxy And Mass Assembly (GAMA) survey 
jPriver et al.l201 lh . lPrescott et alj l l201 lb studied the fraction of red 
satellites around massive galaxies at 0.01 < z < 0.15. They found 
that the fraction of red satellites is higher around massive galax- 
ies than in more isolated regions. Therefore, the satellites would 
be redder in average than the field galaxies in their study, in good 
agreement with our results. 

We also find clear differences between the satellites and their 
central galaxies when we split the sample depending on their Sersic 
indexes. Within the error bars, the satellites of disk-like galaxies 
present a similar evolution of the colors with redshift than their 
central galaxies. In addition, there is no clear differences in their 
mean ages. However, if we focus in the u' - r' color, the satellites 
of spheroid-like galaxies seem to be quite different, since they are 
clearly bluer than their central galaxies and redder than the field 
galaxies. Moreover, the spheroid-like galaxies become significantly 
older than their satellites at low redshift (Fig. [3j. 

The color-magnitude diagrams presented in Fig.[2]suggest that 
at z > 1 the satellites of the spheroid-like massive galaxies could 
be future candidates for mergers involving gas-rich satellites more 
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Figure 3. Difference in stellar age between the central galaxies and their satellites. In the left panel, we plot the case of the whole sample exploring satellites 
down to 1:10 mass ratio (green squares) and down to 1:100 mass ratio (magenta triangles). Results for spheroid-like (red circles) and disk-like (blue stars) 
central galaxies are presented in the right panel. 



than just dry mergers, while at lower redshift we find that satel- 
lites would be more lo cated towards the red s equence. These results 
are in agreement with lNewman et all J2012h . who indicated that a 
significant fraction of mergers at z > 1, even for red hosts, are 
not completely dry (i.e., mergers of massive spheroid-like galaxies 
with disk-like satellites are observed), altho ugh they did not find 
it at lo wer redshifts (0.4 < z < 1). However. lLopez-Saniuan et al.l 
(201 ]]) also found that 35% of their massive early-type galaxies 
(red galaxies) since z ~ 1 present a blue companion, and there- 
fore, they would represent a future mixed merger. Although the 
mean values of the satellites of the massive spheroid-like objects 
seems to be more located in the red sequence at 0.75 < z < 1.10, 
our color-magnitude diagram shows that they are significantly bluer 
than their central galaxies. 

When investigating the ages of the satellites in comparison 
with their central galaxies (Fig. [3j, we find that the satellites have 
similar ages than their hosts at high redshift (z > 1.5), and there- 
fore both massive and satellite galaxies at those re dshifts would 
have formed the bulk of their stars at z > 2 (e.g., Tho mas et al.l 
2005). For galaxies at lower redshifts, we find that satellites are in 
average ~ 1.5 Gyr younger than the m assive galaxies at z = 0. This 
is contrary to the results presented by Pasquali et al. (2010) in the 
nearby universe, who did not find significant differences in the age 
of their satellite and central galaxies with stellar masses as the ones 
considered in this study. However, their classification of satellite 
and central galaxies lies on the identification of the m ost massive 
(or mo st lu minous) group membe r as the central galaxv. lYang et al.l 
( 20051) and ISkibba et al] d201 lh showed that a significant frac- 
tion of central galaxies might actually be a satellite galaxy. Due 
to this potential uncertainty, the results by IPasquali et alj d201 Oh 
could be underestimating the true differences between the prop- 
erties of centrals and satellites. Furthermore, massive galaxies at 
higher redshift (up to z ~ 2) also present n egative color gradients 
(bluer in the outskirts of the galaxies, e .g., Moth & Elston 2002; 
IWu. Shao & etai]|2005l: iMcGrath et alj l20oi ; van Dokkum et al. 



l20ld : lGuoetalJl20TTI ; ICassata et al.ll201ll ; iGargiulo et al1l2012l) . 
This also points out to a picture of younger and/or more metal-poor 
galaxies forming the envelopes of massive galaxies after being ac- 
creted. 

Are the imprints of the high-z satellites found in the outskirts 
of the present-day massive galaxies? Since we find that satellites 
are younger than their centrals at lower redshifts, we should ex- 
pect an age gradient in these massive galaxies. We note, however, 
that due to the small difference at all redshifts between the satel- 
lites and the central galaxies, the age gradient (if an y) would be 
very m ild. Actually, an age gradient was suggested bv lTortora et al.1 
within the errors in their measurements, although other 
authors (e.g., ISanchez-Blazquezet ai] |2007l; iGreene et alj 1201 2|) 
found very shallow or not age gradients. From the well-known 
correlation of stellar mass (or velocity dispersion) with metal- 
licity (e.g.JSanchez-Blazquez et al]|2007l : fkuntschner et alj|2010l : 
Spola or et alJl2010l:lTortora et alj|2010l) 7 our satellites should have 
lower metallicities than the massive galaxies, and therefore, neg- 
ative metallicity gradients would be expected in massive galax- 
ies as it is alread y found (e.g., ISanchez-Blazquez et alj 20071 ; 
IGreene et alj2012l) . 

In addition, since our satellites have lower stellar masses than 
the cen tral galaxies, a more extended star formation history is ex- 
pected dThomas et al.1 1 2005). This would imply lower [a/Fe] val- 
ues for the satellites than for the central galaxies. While massive 
elliptical galaxies in the nearby u niverse have high central val- 
ues of [a/Fe] (~ 0.4 dex; see e.g., ISanchez-Blazquez et al.1 120071 : 
iKuntschner et alJ|2O10t IGreene et al. 2012), less ma s sive g alaxies 
have [o'/Fe]~ 0.2 dex. For example, ISpolaor et al. I d20ld) found 
[a/Fe]~ 0.15 ± 0.11 dex for el liptical galaxies with 3 x 10 10 M o 
(see also lKuntschner et al.|2010l), and late-typ e galaxies could have 
even lower values (e.g., iGanda et ail 120071) . Then, it is reason- 
able to assume that our satellites have similar abundance ratios, 
i.e, [ff/Fe]< 0.2 dex, and we would expect to find slightly nega- 
tive gradients of [ff/Fe] in the nearby massive galaxies. Recently, 
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iGreene et alj j2012h have shown that the stars in outer parts of their 
massive galaxies are metal-poor and [ff/Fe] -enhanced (~ 0.2 dex), 
as it would be expected if our satellite population is eventually ac- 
creted by their central galaxies, suggesting that the outer parts of 
these galaxies are built up via minor merging with a ratio of ~1:10. 
The results presented here also support this scenario and they are a 
benchmark for future detailed studies of the outskirts of local mas- 
sive galaxies. 
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